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ABSTRACT 

Semileptonic B decays are described by the Isgur-Wise form 
factor to the leading order in 1/m; four new functions appear 
in the first order ||] . Values of these functions are crucial for the 
applicability of the whole approach to processes involving c quark. 
We obtain the sum rules for three subleading form factors from 
the QCD sum rules with finite masses by expanding to the first 
order in 1/m. The results respect the pattern of the first 1/m 
corrections established in HQET, and obey the Luke's theorem. 
The numerical estimates show that l/m^ corrections are sizable 
but not catastrophic. 
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1. Introduction 



Recently a significant progress has been achieved in the heavy quark physics 
in the framework of the Heavy Quark Effective Theory (HQET) , see also 
the reviews and references cited therein. The matrix elements of the vector 



current = 07^6 for B — > D^*) decays 



{D*K\B) 



(1) 



to the leading order in 1/m are expressed via the Isgur-Wise form factor: 
^+ = = ^(ch(^, ^_ = 0, where ch(p = vv' , v, v' are 4- velocities of B, 
D^*\ ^(1) = 1 11:01 • First 1/mc corrections involve four new functions 
1/m;, corrections contain no new elements 



1 



1 

1 

mb 
1 



Pi 



1 

rub 

e 

-)- 

rub J 2 



Pi 



(2) 



P2_ 



Pi ~ Pi 

nib 



where the HQET ground state meson's energy s = uib — mt = mn — iric = 
m£)» — rric (these differences are equal up to 1/m corrections). The Luke's 
theorem [5-7] states that pi(l) — 0, ^2(1) — 0. 
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The most interesting physical appUcations of HQET are those to b ^ c 
weak transitions. Therefore the appUcabiUty of the whole approach to pro- 
cesses involving c quark (i. e. the size of 1/mc corrections) is a crucial question 
of HQET. If these corrections are large, HQET has only a purely academic 
interest; if they are modest, it can be applied to real-world processes. 

A nonperturbative method is needed to calculate form factors. QCD sum 
were used for investigation of the form factors (|l|) at finite mb, rric 
HQET sum rule for the Isgur-Wise form factor was considered 
It coincides with the limit m^.c — > oo of the QCD sum rules. 
Results for finite and infinite masses arc compared in p^ . There are two 
alternative ways to obtain sum rules for the subleading form factors pi. One 
can expand the known finite-mass QCD results to the first order in 1/m. 
Alternatively, one can start from the HQET lagrangian and currents in the 
first order in 1/m. The second way gives more insight into the sources of 
the heavy quark symmetry breaking (e. g. the heavy quark chromomagnetic 
moment vertex); the general theorems of HQET (like @ and the Luke's 
theorem) can be traced in the calculation. But the first way is less labour- 
consuming provided that the QCD results are already known, and allows to 
consider easily also higher 1/m corrections. It also gives a strong check of the 
QCD results. Here we use this way. 

The situation is similar in a simpler case of 2-point sum rules. Here the 
QCD (Borel-transformcd) sum rules [15-19] coincide with the leading-order 
HQET sum rules in the limit nib^c oo. The first 1/m correction 

was obtained by expanding the QCD results p^,and also in the framework 
of HQET H. 



2. Sum rules 

We consider three-point correlators 

Kf,{pb,Pc) - J dxbdx,e-'P'^'+'P'^' {TjB{xb)V^{0)j+{x,)) 

= K+{plplt)p^+KApl,pl,t)q^, (3) 
Kf_,^{pb,Pc) = Kv{pl,pl,t)iei_,i,aiiPaqi3, 

where K^^, is similar to with Jd^u instead of jn; jb — q'jsb, jo = 975c, 
jo'iy = qiuC are the currents with the quantum numbers of B, D, D*; p = 
Pb + Pc- q — Pb — Pc, t — q^. We have calculated perturbative spectral 
densities pi{sb,Sc,t) and quark condensates' contributions KKpl,p^,t) up to 
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dimension 6 for the invariant functions Ki{pl,p^,t) using REDUCE |Q; part 
of the resuhs was pubhshed in [ p^ . 

In order to consider the hmit mb,c — > oo, we proceed to the new variables 
pI c — ™6 c + '2'mb,c(^b,c, t = ml + rri^ — 2mbmcch(p. In these variables the 
support of perturbative spectral densities Q is the wedge 

< iUb/iOc < e^. (4) 

At t > spectral densities are singular at the parabola s1 + sj.+t^ — 2sbSc — 
2sbt — 2sct — which touches the boundary of the physical region (^) at 

Lli^ ^ [nib -nice ^) , uJb ^ ^ [-nic + nibe^) . (5) 

The double dispersion representation should be modified above this point . 
In the limit nib^c ~* oo, this area goes to infinity and hence lies outside the 
lowest mesons' duality region essential for the sum rules. 

Making the double Borel transform from the variables uJb.c to Eb^c, we 
obtain QCD sum rules (for finite mb,c) 

^^C±(ch^)e-(-+-)/(-) = 
nib rric y/msmD 

I J [{mB ± mD)p+ + (me T mD)p-] e-^''-+'^^^'^^^UuJbduJc 

+ AE^B[{mB±mD)Kl^{mBTmD)K'L\^, (6) 

^/D.mz,.Cy(ch¥.)e-(--+--)/(2^^) 

7 jy" pye-('^''+'^^)/(2£)^^^^^^ ^ AE-'BK. 
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where we have taken Eb = Ec = 2E; m\ = m^ + 2mt,eB, and similarly for 
D^*K The integrals are calculated over the part of the wedge ^) dual to 
the lowest mesons in both channels, i. e. the wedge minus the higher states' 
continuum region. It is convenient to introduce the variables w, rj instead 
of u)b,c = uj [l ± ijth^ ^) so that the physical region (|4[) is — 1 < ry < 1, and 
to measure quantities with the dimension of energy (such as E) in units of 

k, = ^ \ {qq)\ (mo — ikEo). Expanding to the first order in l/mb^c, we 
obtain 



(m6mc)3/2| (qq) \ 
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A{^,E)+[— + —]kB,{^,E), 
\mc mh' 



(m{,mc)'''/2| (qq) \ 
1 1 



-A{ip,E) + kBi{(p,E) 



(7) 



ml''^rnj^\ {qq) \ 
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16 ch^ ^ 



3 E^ 277r£;3 
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8ch^f 



(3 - r]^)Lo^e~'^/^dr^dLu 



Similar expansion of 2-point sum rules gives 

/n"^n . /P ^ 2kBiiO,E) 



277r£;2 ' 
2El as{Ac\Yip + l) 



27ttE^ 



■ml \ {qq) \ 

"T-cl (gq) I 



rric 
2kB2{Q,E) 



(9) 



(the sum rule for fs is similar to the one for fo). From (0) and (|) we obtain 
the sum rules for ^(chi^a) p^,p3[ and for /Oi,2.4(ch (ys) 



, , A(w,E) , Bi(ip,E) BAQ,E) , , 

g(chv?)^ / , pi(chy)/fc^ / - An J . (10) 



^(0,S) 



^(^,S) ^(0,S)' 
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B2{ip,E) B2{0,E) Bi{^,E) 

The results of 1/m expansion of the QCD sum rules obey the structure (Q) 
and the Luke's theorem. It is easy to explain why pi(l) = 0. At finite 
nib = q the 3-point correlator is related to the 2-point one by 
the Ward identity K+{p^,p^) = ^21p!) ji^^ and hence the QCD sum rule 
reproduces the exact result = 0) = 1. At infinite m^.c, the Ward identity 
is K{lo, lo) ~ '^^^^ , and hence the HQET sum rule reproduces the exact result 
^(1) — 1. Then we obtain from (||) pi(l) = 0. Therefore it is not accidental 
that the 1/rn correction to the sum rule for (|^) is described by the same 
function Bi as to the sum rule for ^+ but at 1^9 = 0. The similar fact for 
the sum rule for Jd- and l/rric correction to (,v (0) looks like a miracle in 
the framework of QCD sum rules, but it leads to the Luke's theorem for p2 
in (|l0|). The dependence between l/rrif, correction to i^y and 1/m corrections 
to ^± in (j^) is one more miracle. From (j^) we obtain (compare with [2^ ) 

fo = ^{l + ^), fo^^^(l + ^), (11) 



where the coefficients are given by the sum rules 

^e-."= = ^(O.B), (12) 



3. Results 

At the standard values of the condensates, k = 260-280MeV, Eq = ^ = 
0.85. The 2-point sum rule ( [l^ ) for was analyzed in j2^. The continuum 
threshold is o^o = 3 and the meson energy is e = 1.65 in the Borel parameter 
plato E = 1.7-2.5 (all in k units). Now we know that the perturbative 
correction to these results is large But this correction is unknown in 

the 3-point sum rule. We may hope that it will be partially compensated in 
the ratio of 3-point to 2-point sum rules. For consistency, we use the 2-point 
results without perturbative corrections as an input for the 3-point sum rules 
analysis. In any case, even a rough estimate of pi is currently valuable. 

In the case of 3-point sum rules, e. g. for ^{ch ip) (|o|), the higher states' 
continuum spectral density is modelled by the perturbative spectral density 
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Figure 1: Continuum models 



with some smooth curve as a continuum threshold (Fig. At </? — > the 
wedge becomes infinitesimally narrow, and the continuum starts at coq for 
consistency with the 2-point sum rules. The wedge area is f, therefore the 
spectral density must be ~ 1/ip in order to yield the perturbative contribution 
~ 1 (of course, the explicit calculations confirm it). The simplest continuum 
model is that with the straight line threshold. If one chooses any smooth curve 
instead, the difference area is ^ ip^ (Fig. [l]), and the variation of ^(chc/?) is 
^ (fi'^. It influences the slope of £,{ch ip) at 1; this freedom is analogous to the 
freedom of choosing the continuum threshold in the 2-point sum rule (|l^ ) for 
But if one chooses a continuum threshold with a cusp in the physical 
region (as was done in fl^), the difference area is ~ (p'^ (Fig. [|), and the 
variation of (,{ch(p) is ~ ip. This gives an infinite slope of £,{ch ip) at 1, 
what contradicts to its general analytical properties. Besides that, there are 
no physical reasons for the continuum threshold to be non-smooth. Here 
we restrict ourselves to the simplest continuum model with the straight-line 
threshold; dependence of f(chiyj) on the threshold curvature was investigated 
in(l|. 

The sum rule for ^ (ch ip) is well known [T^ , p^ . It contains two main 
terms: the perturbative contribution and the quark condensate one. The 
quark condensate contribution is constant (up to nonlocality effects discussed 
in Ip^). The perturbative contribution contains the light quark propagator 
from Xc to Xb, and falls as 1/ {xb — Xc)'^. The Euclidean distances from to Xb 
and from to 1/(2_E); the angle between these lines is (p. Hence the 

distance from Xb to Xc is ch ^ /E. The perturbative contribution is suppressed 
as l/cli^j compared to the 2-point (straight-line) case. This leads to the 
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decreasing of ^(ch(^) with increasing Lp. Graphs for ^(ch(^) at several Borel 
parameters are shown at Fig. 0; the remarkable stability is seen. 




0.2 - 



1.25 1.5 1.75 2 2.25 2.5 2.75 3 



Figure 2: Isgur-Wise form factor ^(ch(p) at Borel parameter values E = 1.7, 
2, and 2.5 

Now we present the new results for the 1/m correction form factors 
Pi.2.4(ch (^). It is seen from the sum rules ( p^ that their scale is set by 
/ 2 _ \i/3 

the unit k = [^\ (qq) \ j — 260-280MeV. The dimensionless functions in 
the sum rules don't contain large numerical factors. The numerical analysis 
of the sum rules confirms this conclusion; the results are shown at Fig. |. The 
curves pi.2(ch(/3) are pinned at the origin by the Luke's theorem; at cht/? > 1 
they grow and reach the values of order k. On the other hand, pi^chip) is 
small and nearly constant. Variation of the results with the Borel parameter 
E allows to estimate their accuracy. 

In this work we have investigated the sum rules for pi^2,4(ch (^). It is 
sufhcient for the decay B D. Form factors of the decay B D* also 
contain p^. In order to obtain sum rules for it, it is necessary to consider 
an axial correlator in addition to (^. Such an analysis will be presented 
elsewhere. 

When this work was completed and presented at the seminar at SLAG, 
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2r 




Figure 3: 1/m correction form factors pi{c]iip) (solid curves), p2{chip) 
(dashed curves), and pi{c]iip) (dashed-dotted curves) at Borel parameter val- 
ues E = 1.7, 2, and 2.5 

M. Neubert informed us about his preprint | p6[ . In this work the sum rules 
for the subleading B D^*^ form factors were obtained in the framework of 
HQET. The results are similar to our ones. We are grateful to M. Neubert 
for giving us the preprint ||2^ and for the useful discussion. 

A. Correlators in QCD 

Here we present the spectral densities and quark condensates' contributions 
up to dimension 6 for the correlators (^ and the correlator similar to K^j.^, 
but with the axial current c'jfj.^sb (it has the structure KAgp.^ + K++p^pu + 
+ K 

— \-QiiPv + K — qiiqv)- The results for all correlators have been 
produced by a single REDUCE program in which only few lines with 
the 7-matrix structures of currents varied, this allows to avoid bugs in the 
programs for separate channels. To minimize the probability of printing er- 
rors, a single REDUCE source was used both for algebraic checking and for 
production of I^TJt;X source of equations in this appendix (using the package 
RLFI by R. Liska from the REDUCE library H). 
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The spectral densities are 
N 



P+ 
P- 

Pv 

PA 
P++ 

P+- 



87r2A3/2 

N 



t{2xbXc + {rub + mc)a+) + (rnt - mc)a+h 



- {{mb - rricY - t) {mb - mc)a+ + 2xbXc{xb - Xc) 



87r2A3/2 

+ a_((3mc — mb)xb + (3mfc — mc)xc) 
N 



87r2A3/2 
N 



87r2A3/2 

NrUb 



P-+ = 



{{mb - mcf -t)a+ + {xb - Xc)a 

Aa+ + 2mb{xbXct + {xb - Xc){mlxb - mlxc)) 

g^2^5/2 I'^^i^b + Xc) + A{6xbXc - {ml - ml){xb - Xc)) 
+ 6t{xb + Xc){2xbXc + m^Xb + mixc) — Qxbxjp' 

- A.{(^Xb - Xc)a+ - {mb - mc){{ml + ml)xb - 2m; 

- 6mb{xhXct + {xb - Xc)imlxb - mlxc))b , 
8^2A5/2 [-Ateb(mb + me) 

+ A({xb - Xc)a- + {mb + mc)((mj + ml)xb - 2mbXc)^ 

- 6mb{xbXct + {xb - Xc){mlxb - mlxc))b 



I'b-^c, 



Nmb 



At{xb — Xc 



87r2A5/2 

+ A{-2xc{2xb - Xc) - i'Sml + ml)xb + (3m^ + ml)xc) 

- 6t{xb - Xc)imlxb - mlxc) 

- 6{xbixb - Xc) + {ml + ml)xb - 2mlxc) 

{xc{xc - Xb) + {ml + ml)xc - 2mlxb) , 



where A 



„2 I „2 I +2 



s'^^+V-2sbSc-2sbt-2sct, Xb,c = Sb,c-m.l^, a± = mcXb±mbX. 



b = Xb — Xc + ml — ml- 

The results for g(75,7i/)c ^ 7^.75)0, 6(7m,7a«75)c ^ &(7/:i75, 7a«)c can 
be easily obtained by 

= -Kt{mc -mc). 
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There is an interesting check of these formulae. If we multiply a correlator 
with g7^c by Pci, and use the identity Sq{k)'yuSc{Pc + k)pcv = mcSq{k)Scipc + 
k) + Sq{k) — Sc{pc + k), we obtain the corresponding correlator with qc plus 
terms having no double discontinuity. Therefore 

PA + [sb + isc -t)p++ + {sf, - Sc -t)p+^ = -2p+{mc ^ -rric), 
-pA + (sb + 3sc - t)p-+ + {sb - Sc- t)p — = -2p-{mc ^ -rric). 

The quark condensates' contributions arc 
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i{Q9GlJ"(Ti^^q) / {qq), m? 
factorization approximation (J^ 
A'^ = 3 is the number of colours). 
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